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Abstract 

CN ' A new type of an integrable mapping is presented. This map is equipped 

| with fractional difference and possesses an exact solution, which can be re- 

garded as a discrete analogue of the Mittag-Lefner function. 

Hi 

m ■ 1 Introduction 

Logistic map |5|, U 

GO . 

C; — = au n (l - u n+x ) (a > 0) (1.1) 

e 

is an integrable discretization of the well-known logistic equation. 

^-u = au(l-u) (a>0) (1.2) 

Its origin was a population model in ecology. In addition, it has been reported U 
that the logistic map is used in many other fields, for example, agriculture, life 

(N 
O 



X 



sciences and engineering. 

Through dependent variable transformation, 



C ■ l + ae 

u n = (l.dj 

> . aeg n + 1 + ae 



the logistic map ( |1.1|) is linearized as 

A_ n # n = -ag n , (1.4) 
where A_ n is a backward difference operator defined by 

A_„ = e-^l - E- 1 ), E~ x f n = / n _x. (1.5) 
Since equation ( |1.4j ) possesses a solution, 

g n = g (l + ae)- n (1.6) 
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a solution to the logistic map ( |1.1|) is given by 



u + (1 - u )(l + ae)- n V ; 

The main purpose of this paper is to provide a new type of an integrable mapping, 
which can be regarded as an extension of the logistic map and is equipped with 
fractional difference. This map, which we call the fractional logistic map here, 
possesses an exact solution and have another parameter p which corresponds to an 
order of difference. In section 2, we introduce one definition of fractional difference, 
which is a slight modification of Hirota's fractional difference operator. We also find 
an eigenfunction of this operator. Section 3 is the main consequence of this paper 
and presents the fractional logistic map. We also show its time evolution through 
numerical experiment. 



2 Fractional difference 

We here give a definition of fractional difference operator and its eigenfunction. 
Before going to its definition, let us introduce fundamental functions M(a; n) defined 
by 

1 „ i Tin + a — 1) n i I n + a — 2 \ 
M{a-n) = —e a - 1 K = e a ~ l o>0,nGZ, 2.1 

1 (a) L(n) \ n — 1 / 

where e is an interval length and I ] (a £ R, n G Z) stands for a binomial 

y n 

coefficient defined by 

' a(a- !)■■■(« -ra + 1) r(a+l) 



n\ r(a-n+l)r(n+l) ^ > °^ 

1 (n = 0) 

(n < 0) 

This function satisfies the following lemma. 
Lemma 2.1 The following relations hold. 

M(a; 0) = (a > 1) (2.2) 
A_ n M(a + l;n) = e _1 (M(a + 1; n) - M(a + l;n- 1)) = M(a;n) (a > 0) (2.3) 
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Proof: Equation ( |2.2|) is obvious owing to the definition of binomial coefficient. 
Equation ( p. 3D is proved by using the relation T(x + 1) = xT(x) as follows. 



e~ 1 {M{a + 1; n) - M(a + 1; n - 1)) 



-a-1 



[n + a- l)T(n + a- 1) (n - l)r(n + a — f 



r(a + 1) v r(n 
ae^ 1 r(n + a - 1) 



rfn 



r(a + i) 



r(n) 



M(a;n). 



Next we go to the definition of fractional difference. Hirota Q took the first n 
terms of Taylor series of A" n = e~ a (l — E^ 1 ) 01 and gave the following definition. 

Definition 1 Let ael Then difference operator of order a is defined by 

-l)% n _3 a ^ 1,2, • • • 



n-l 



a 



A?„«„ = •! 



3=0 \ J 

m 



(2.4) 



-m 



77? 



^-l)% n _j a = m G Z >0 



i=o \ J 

It should be noted that Diaz, Osier JJ] gave another definition of fractional difference, 



n-j 



(2.5) 



3=0 



We here adopt another difference oparator A"__ n by modifying Hirota's operator. 

Definition 2 Let a G R and m fre an integer such that m — 1 < a < m. VFe define 
difference operator of order a, A"_ n; fry 



n-l 
3=0 



We define a new function, 



a — m 



•l)^AVfcU=n- 



(2.6) 



F a (A,n) = ^ AW(aj + l;n). = ^A% aj ' : 



3=0 



3=0 



T(n + aj) 
T(aj + l)r(n) 



(2.7) 
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Remark 1 Putting a = 1 in the above definition, we have 

T(n+j] 



CO 



T(j + l)r(n 

OO 

= t» 



j=0 

n + j — 1 



J 



5>A £ y( ^. )=(l-A£)-». (2. 



j=0 



Remark 2 In t/je limit of e — > 0, n — > oo wift t = ne fixed, the function F a (X,n) 
converges to 

E a (x) is the well-known Mittag-Leffler function and and its asymptotic behavior is 
studied in detail in ffij. 



Theorem 2.1 If a > 0, the function F a (X,n) is an eigen-function of the fractional 
difference operator ( |2.6| ). That is, 

A a F a (X,n) = XF a (X,n) (2.10) 



Proof of Theorem |2.1| : Let m be an integer such that m — 1 < a < m. Then we 
have 



Al_ n F a (X, n) = A*_ n 1 + 53 X j M(aj + 1; 



A: i _„53AW(aj + l; 



oc 



53 X j A a _- n m M(aj + 1 - m; n) (2.11) 
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Each summand in the above equation is given by 

A a _- n m M(aj + l-m;n) = ^[ _ j ( _1 ) fcM ( a i + l-m;n-k) 



fc=0 

n-1 



El a — m \ . k I aj — m + n — k — 1 
l u I ^ l ^ _ h _ 1 



fc=0 
n-1 



A; y \ n — A; — 1 

■1)- f - OJ - 1 ) = f ° J - ' - 1 | = M(aj - . + 1; „). (2.12) 



Therefore, substitution of eq. (|2.12|) into eq. (|2.11|) gives 



A* _ n F a (A, n) = J2 XM{aj - a + 1; n) = A i+1 M(aj + 1; n) = AF a (A, n) 

3=1 j=0 

(2.13) 

which completes the proof. ■ 



3 Fractional Logistic map 

This section presents a fractionl integrable mapping, which we call fractional logistic 
map here. We start with a linear equation 

K,- n 9n = ~ag n (0 < p < 1, a > 0) (3.1) 
instead of linear eq. ( |1.4|) . The above equation is rewritten as 



g n = (l + ae?)- 1 j<?n-i - £ y j 1 J (-^(^-i " ft-i-i)} (3.2) 

Through the similar dependent variable transformation as logistic map, 

u n = 3.3) 

aePg n + 1 + aeP v 7 

we finally obtain a fractional logistic map 



Due to the Theorem p7T| , eq. (|3.1|) has a solution, 



9n = g F p (-a;n). 

Therefore, a solution to the fractional logistic map ( |3.4| ) is given by 

uo 

u + (1 - u )F p (-a;n) 

Putting p — 1 in eq. ( |3.4|) , we have 

(1 + CLE^Un-X 



(3.5) 



(3.6) 



u. 



1 + (1 + ae)u n ^i 



(3.7) 



which recovers the logistic map. The following figure illustrates time evolutions of 
fractional logistic map with order parameter p = n/4(n = 1,2,3,4). We have put 
Mo = 0.1, a = 1.0 and e = 0.1. 




Figure 1: Time evolutions of fractional logistic map 



Considering the fact that the Mittag-Leffler function has an asymptotic behav- 



ior 



1 \-7i4-— an 



■i T(l — an) 



0(t 



-aN\ 



t — > oo, A < 0, 



{3.1 
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we can observe that u n converges to 1 at the order of 0(l/n p ) if < p < 1. The 
following table illustrates a numerical result in which we apply the p-algorithm [l] 

(k + n) p - n p 



Pk+1 



P-i = o, Pi 



Pl +1 



Pi 



to the sequence {u n } in the case p = 1/4. This table shows that u n converges to the 
value near 1.0 at the order 0(l/n 1/4 ) as n tends to +oo. 



Table 1: The p-algorithm applied to the sequence {u n } in the case p = 1/4 





Pi 


Pi 


Pi 


P7 




P21 


1 


0.1 


0.19745 


0.84609 


0.84288 




1.00145 


2 


0.14791 


0.23921 


0.86006 


1.30904 




0.99915 


3 


0.16019 


0.27802 


0.99118 


1.21129 




1.00131 


4 


0.16764 


0.31406 


1.07379 


1.17824 




1.00139 


5 


0.17313 


0.34778 


1.11866 


1.16134 




1.00112 


6 


0.17752 


0.37941 


1.14064 


1.15284 




1.00120 


7 


0.18121 


0.40908 


1.14968 


1.15045 




1.00120 


8 


0.18442 


0.43692 


1.15161 


1.15375 




1.00121 


9 


0.18726 


0.46304 


1.14970 


1.16376 




1.00122 
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